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Dark matter decays or annihilations that produce line-like spectra may be smoking-gun signals. 
However, even such distinctive signatures can be mimicked by astrophysical or instrumental causes. 

We show that velocity spectroscopy—the measurement of energy shifts induced by relative motion 
of source and observer—can separate these three causes with minimal theoretical uncertainties. The 
principal obstacle has been energy resolution, but upcoming experiments will reach the required 0.1% 
level. As an example, we show that the imminent Astro-H mission can use Milky Way observations 
to separate possible causes of the 3.5-keV line. We discuss other applications. 


Introduction: What is the dark matter? Identifica¬ 
tion depends upon more than just observation of its 
bulk gravitational effects; distinct particle signatures are 
needed. Backgrounds make it difficult to pick out these 
signals, which are constrained to be faint. Among possi¬ 
ble decay or annihilation signals, those with sharp spec¬ 
tral features, such as a line, are especially valuable. 

Given that the stakes and difficulties are so profound, 
even such a “smoking-gun” signal may not be conclusive. 
A line could have other causes: astrophysical (baryonic) 
emission or detector backgrounds (or response effects). 
For example, the cause of the recently discovered 3.5-keV 
line is disputed [1-8]. This problem is more general [9- 
18] and will surely arise again. We need better evidence 
than just a smoking gun—we need to see it in motion. 

Premise and Motivation: We propose a general 
method for distinguishing the possible causes of a sharp 
spectral feature. Gonsider a line of unknown cause—dark 
matter (DM), astrophysical, or detector—observed in the 
Milky Way (MW). Relative motion between source and 
observer leads to distinctive energy shifts as a function 
of line of sight (LOS) direction. Figure 1 illustrates this 
schematically. Because typical Galactic virial velocities 
are ~ 10“^c, the Doppler shifts are only ~ 0.1%. 

A potential target for velocity spectroscopy is the 3.5- 
keV line recently observed in MW, M31, and galaxy clus¬ 
ter spectra [1, 2, 4]. The line energy and flux can nat¬ 
urally be explained by sterile neutrino DM [19-28] (or 
other candidates [29-37]). However, the significance of 
the line is disputed [3, 5, 6], and it has been argued that 
it can be explained by astrophysical emission [7, 8]. 

With present detectors, velocity spectroscopy of this 
line is impossible. Excitingly, the Soft X-Ray Spectrom¬ 
eter (SXS) on Astro-H (launch date early 2016) has a goal 
energy resolution of ctah = 1-7eV (4 eV FWHM) [38, 39], 
which is at the required 0.1% scale. We show that if 
this goal resolution is achieved, Astro-H can identify the 




FIG. 1. Top: How DM, astrophysical, and detector lines 
shift with Galactic longitude is starkly different. Bottom: 
For DM signals at positive longitude, our motion through 
the non-rotating DM halo yields a negative LOS velocity and 
thus a blue shift. In contrast, for astrophysical lines (e.g., 
from gas), co-rotation in the disk leads to a positive LOS 
velocity and thus a red shift. These signs reverse at negative 
longitude. Detector lines have zero shift. 


cause of the 3.5-keV line. We also discuss prospects if 
the performance is worse. 

We emphasize that the applicability of DM velocity 
spectroscopy is much more general. The purpose of this 
paper is to introduce a new concept to increase the power 
of DM searches and to spur innovation in detector design. 
We conclude by discussing several generalizations. 
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Usual DM Decay Signal: The differential intensity 
(flux per solid angle) from DM with mass and lifetime 
T = 1/r, decaying within the MW, is 


dE 
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where Rq ~ 8 kpc and Pq ~ 0.4 GeV cm“^ [40-42] are 
the distance to the Galactic center (GG) and local DM 
density. (We neglect the cosmologically broadened extra- 
galactic signal, which contributes negligibly in Astro-H’s 
narrow energy bins.) i/('0) is the dimensionless, astro- 
physical J-factor defined by the LOS integral 


J{tp )=^— [ ds p^{r[s,ilj]), (2) 

Rq Pq J 


where ip is the angle relative to the GG and is related 
to Galactic longitude and latitude via cos ip = cos I cos b. 
dN{E)/dE is the photon spectrum. 

The above treatment assumes that the astrophysical 
term, J{ip), and the photon spectrum, dN{E)/dE, are 
separable. However, for detectors with energy resolution 
< 0.1%, this approximation is not valid because relative 
velocities between source and observer, and therefore the 
spectral shape, vary along the LOS. 

Modified DM Spectrum: We first account for how 
the signal is broadened by DM velocity dispersion and 
second for how it is shifted due to bulk relative motion. 

We take the DM halo of the MW to be spherically 
symmetric, in steady state, and to have no appreciable 
rotation. The last is expected from angular momentum 
conservation, as the baryons from the proto-halo have 
collapsed significantly, while the DM has not; this is con¬ 
firmed by simulations [43, 44]. Thus, (v^) = 0. 

DM particles do have non-zero velocity dispersion, 
determined by the total gravitational potential of the 
halo [45, 46]. Assuming an isotropic velocity distribution 
io'v,r = cTy,!!) = SO the total dispersion is VScfy^r), 

the radial velocity dispersion of DM is [45] 
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where Mtot(r) is the total mass within a radius r. Typical 
values at r ~ few kpc are (Ty^y ~ 125kms“^. 

To calculate ay^r{r), we adopt the mass model of 
Ref. [47], which fits a contracted DM and three- 
component baryon mass profile to MW rotation curve 
data; for more details see Supplemental Materials. The 
choice of mass model is not critical; kinematic results 
from other models agree within 0(10%) [40, 48]. 

The spectrum from a point along the LOS is the con¬ 
volution of the intrinsic spectrum with the DM velocity 
distribution at that point. We assume a Maxwellian ve¬ 
locity distribution throughout the halo, which, at each 


point, yields a Gaussian distribution of the LOS velocity 
component. The modified spectrum from each point is 


dN{E, r[s. Ip]) 
dE 
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where G{E\ ue) is a Gaussian of width ge = {E/c)ay-^y^^. 
Based upon observations of the LOS velocity distribution 
of MW halo stars reported in [49], we take Gy-^^^{r) ~ 
(^y,r{r) which implies ge = {E/c) Gy^r{'r[s,ip]). 

The line shift follows from the LOS velocity, Ulos = 
((^A) ~ do) ■ Tlosi where positive Ulos indicates reced¬ 
ing motion. For Ulos ^ c, the resultant energy shift is 
6Emw/E = —Vi^os/c. 

The Sun follows a roughly circular orbit about the GC 
in the direction toward positive Galactic longitude at a 
speed Vq ~ 220kms“^ [50]. (Recent work suggests Vq > 
240kms“^ [51, 52], which would strengthen our results.) 
The spectrum is therefore shifted by SEMw{hb)/E = 
-|-(uo/c) sin I cos b, which changes sign with 1. We neglect 
the solar peculiar velocity as well as Earth and satellite 
motions, all of which are < 10kms“^ [53-55]. 

The final expression for the modified spectrum, includ¬ 
ing broadening and shifts, is therefore 


d^ _ 1 

dE RqPq 




(5) 

so that Eq. (1) is altered by J (ip) dN {E) / dE -Q 
dj(ip, E)/dE. The observed signal, which is the con¬ 
volution of dJ/dE with the detector response, is nearly 
Gaussian and has an effective width GeS- 


Modified Astrophysical Spectrum: The details are 
slightly different for astrophysical lines. 

The widths of astrophysical lines are primarily deter¬ 
mined by the mass of the emitting atom and by the 
gas temperature; turbulent broadening is negligible [56]. 
Eor potassium at T = 2keV, the intrinsic line width is 
CTgas — 0.8 eV, comparable to Astro-H’s goal resolution, 
ctah — l.TeV. The intrinsic width is weakly sensitive to 
the gas temperature and mass (oc y/Tjrn); any reason¬ 
able values of T and m give similar results. 

For the shift of an astrophysical signal, we must ac¬ 
count for co-rotation within the MW disc. (While there 
is a non-rotating, gaseous halo at the outskirts of the 
MW, it is not hot enough to produce significant emis¬ 
sion at 3.5 keV [5, 57, 58]). For simplicity, we assume all 
baryons follow circular orbits about the GC with speed 
VciTcir) = \/GMtot{r)/r. With this circular speed and 
the hot gas distribution of Ref. [59], we compute the 
spectral shift by integrating the signal along the LOS 
with the contribution from each point weighted by the 
gas density. We call this fiducial model G2. 

Because the spatial and speed distributions of MW X- 
ray gas are uncertain, we compare to models in Ref. [60] 
with smaller and larger line shifts. G1 is based on the dis¬ 
tribution of free e~ [61] and the MW rotation curve [62]. 
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FIG. 2. Comparison of received spectra for DM and gas (G2). 
The emitted spectra are taken to have equal flux and to be 
centered at 3.5 keV before velocity effects. The line profiles 
include velocity dispersion and shift effects, as well as the 
energy resolution of Astro-H. Vertical bands indicate the 1-a 
centroid uncertainties after 2-Ms observations. For contrast, 
the brown line in the figure and inset shows the same signal 
if Astro-H had the energy resolution of XMM. 

G3 is based on the observed distribution of ^®A1 gamma 
rays [60]. G1 and G2 are in good agreement with MW 
HI and CO data [63, 64]. Peak LOS velocities for Gl, 
G2, and G3 are ~ 50,75, and 250kms“^. 

Line Flux Detection: One prerequisite to detecting a 
spectral shift is that the number of signal events be non¬ 
zero. Another is that the background fluctuations be 
small in comparison. Though Astro-H has a small field 
of view (FOV), its excellent energy resolution strongly 
suppresses backgrounds for a line signal, so that even a 
small number of signal events can be significant. 

Viewing directions I ~ 10° — 40° have advantages. 
First, the balance between decreasing signal flux and in¬ 
creasing energy shift at large I is optimized. Second, 
theoretical uncertainties are minimized, as the DM den¬ 
sity profile at r > few kpc is fixed by rotation curve 
data. Third, continuum astrophysical backgrounds are 
reduced; we reduce these further by going slightly off the 
Galactic plane, which minimally affects the DM signal. 

The expected signal intensity is calculated from 
Eq. (1). For our DM example, this is 

I{ijj) = 1.2 X 10“® cm“^ s“^ arcmin“^ (6) 

sin^ 26 > \ / Wx W / J{ip) \ 

7 X 10-11 J V7keVy \J{1 = 20°, [fo] = 5°)/ ’ 

where we have integrated over energy in the line profile. 



FIG. 3. LOS velocity for DM and various gas models (the 
realistic version of Fig. 1). Uncertainties are computed as¬ 
suming 2-Ms Astro-H exposures on each point. 

calculated J{1 = 20°, \b\ = 5°) = 7.5 using Ref. [47], and 
taken the DM parameters from Ref. [Ij. For Astro-H, 
^2 fov = Oarcmin^ and Aegf = 200 cm^ [38, 39], so the 
expected number of events is 

Ns{^) ~ 43 (^j^i ^20°, \b\ =5°)) (^) ■ 

This assumed exposure is large, but appropriate to the 
stakes (a potential discovery of DM) and the difficulties 
(the total exposure of XMM, Ghandra, and Suzaku used 
in the 3.5-keV analyses is > 40 Ms [1, 2, 4, 65, 66]). Fur¬ 
thermore, due to Astro-H’s excellent energy resolution, 
all pointings in a substantial fraction of the sky will help 
test the 3.5-keV line. 

For continuum backgrounds, we consider only the con¬ 
tribution over the narrow energy range ±2CTeff centered 
at 3.5 keV. (We do not need to include the tails of nearby 
astrophysical lines, as they will be well-resolved, unlike 
in XMM.) One component of the background is due 
to the isotropic cosmic X-ray background (CXB) [67- 
69]. We conservatively adopt the total GXB flux (un¬ 
resolved -I- resolved sources) E d^cxs/dE = 9.2 x 
10-^(E/keV)-°'^ cm”^ s”^ arcmin”^ [69]. Another back¬ 
ground, due to hot gas in the MW, varies strongly with 
direction [70]. Finally, there are detector backgrounds 
due to intrinsic and induced radioactivities as well as 
cosmic-ray interactions; their intensity is expected to be 
comparable to that of the CXB [71]. For ■0(1 = 20°, [fe] = 
5°), backgrounds contribute TVf, ~ 5.2 -|- 5.4 -|- 5.4 = 16 
events per 2 Ms within the ±2(Teff — ±4.8 eV band cen¬ 
tered at 3.5 keV, compared to ~ 41. 
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We estimate the detection significance by the Poisson 
probability P{in > 57|/r = 16), which corresponds to a 
one-sided Gaussian probability > la. 

Line Shift Detection: Detecting a line shift depends 
on how well the centroid of the line profile is determined. 
Backgrounds decrease the precision, but, as above, the 
energy resolution of Astro-H plays a critical role. 

When backgrounds are absent, the uncertainty on the 
centroid is creff/-\/iVs- When they are present, the un¬ 
certainty becomes 5E = C{R) aee/\/N^, where C(R) is 
a correction factor and R is defined by the background 
to signal ratio. We calculate the optimal C(R) using the 
Cramer-Rao theorem [72-74]. For '0(1 = 20°,|6| = 5°), 
C{R) ~ 1.6, so that the uncertainty in the LOS velocity 
is — 50kms“^. 

Figure 2 shows the line profiles at 0(1 = 20°, j&j = 5°) 
for a 3.5-keV emission line, due either to DM or gas. 
(A detector line would have zero shift). These profiles 
show how the energy spectra are shifted due to relative 
motion as well as broadened due to intrinsic dispersion 
and detector resolution. We show the uncertainties on 
the centroids, which are separated from each other and 
from zero in a 2-Ms exposure. With the energy resolution 
of XMM [75] (crxMM 47 eV vs. ctah 1.7eV), the 
profiles are indistinguishable. 

Figure 3 shows how the expected shifts vary with 
Galactic longitude, along with their uncertainties, assum¬ 
ing 2-Ms observations for each point. We show the DM 
signal uncertainties; for an astrophysical line of the same 
flux, the uncertainties are comparable because the effec¬ 
tive widths are comparable (a^g^ ~ 160kms“^, — 

200kms“^); see Fig. 2. For a detector line with zero in¬ 
trinsic width, the effective width is ~ 150kms“^, 
approximately a factor of less than 

For each point in Fig. 3, it is easy to assess the prob¬ 
ability that the expected DM signal could fluctuate to 
match that expected for an astrophysical or detector line, 
i.e., that a true DM signal could remain hidden. With 
two observations, at I = ±20°, this scenario can be ruled 
out, relative to G2, at ~ 3.6a. This establishes that this 
technique has interesting sensitivity. Once there is data, 
one can assess the probability that an astrophysical or de¬ 
tector line could mimic a DM signal (for the same flux, 

Agas ^ Adet ^ ADM / /o'! 

“«LOS “ "'I’LOS “ 

If the energy resolution is worse than the design goal, 
e.g., tJAH — 2.1, 2.5, or 3 eV, then the line shift signifi¬ 
cance is ~ 3.0, 2.4, or 1.9cr (the line flux significance is 
always > 5a). This could be improved as -s/t with more 
exposure (including non-dedicated pointings). We have 
not included the systematic uncertainty due to detector 
gain calibration, for which the goal is 0.4 eV [71]. This 
can be mitigated by comparing the energies of nearby 
astrophysical lines, especially at opposite longitudes. 

Related Searches: Astro-H may be able to resolve the 
intrinsic width of a MW DM line. This would provide the 


first information on the large-scale DM velocity distribu¬ 
tion, which is sensitive to DM particle properties [76] and 
to the presence of substructure [77, 78] (see Suppl. Mat.). 

The 3.5-keV line has been detected in M31. Due to the 
relative motion between the Sun and M31, DM or astro- 
physical lines from the center of M31 will have LOS shifts 
of ~ — 300kms“^ [79]. We estimate that this blue shift 
could be detected with > 5cr significance, making this 
an attractive way to test detector causes. Due to M31’s 
rotation, astrophysical lines are separated from DM lines 
by ±200kms“^ around ±1°, but, because the statisti¬ 
cal uncertainties are large, they cannot be cleanly distin¬ 
guished in 2 Ms; see Suppl. Mat. and Refs. [80-86]. The 
LMC [87] may also be an attractive target. 

More speculatively, it may be possible to see the line in 
the extragalactic DM signal, if more astrophysical sources 
in the GXB are resolved, e.g., with eRosita [88, 89]. Fur¬ 
thermore, because we move at ~ 400kms“^ with respect 
to the CMB, it may be possible to detect a dipole signa¬ 
ture in DM line signal. Far-future observations may even 
detect a forest of sources in each LOS spectrum. 

Conclusions: Even for a supposedly smoking-gun sig¬ 
nal, such as a line, it may be difficult to distinguish be¬ 
tween DM, astrophysical, or detector causes. We have 
shown that detectors with energy resolution < 0.1% can 
break this degeneracy using velocity spectroscopy, which 
has minimal theoretical uncertainties. We emphasize 
that our main goal is to point out this new and robust 
method for testing DM signals, which can be applied to 
any sharp feature, such as an edge or box [90, 91]. 

To demonstrate the potential of this technique, we have 
shown that Astro-H will be able to test the origin of the 
3.5-keV line. In the future, other lines may be discovered. 
For lines at higher energy, the relative energy resolution 
of Astro-H improves. This unprecedented resolution will 
allow Astro-H to dramatically improve on existing ster¬ 
ile neutrino limits [92-110]. We encourage a dedicated 
study by the Astro-H Gollaboration, once post-launch 
parameters are known, to give definitive answers on DM 
sensitivity over their full energy range. 

We are encouraged by the expected 0.1% resolution of 
Astro-H in the range 0.3—12keV, and the demonstrated 
0.1% resolution of INTEGRAL-SPI in the range 20keV 
to 8MeV (including velocity spectroscopy of the 1.809- 
MeV line from ^®A1 [60, 111, 112]). Excitingly, the pro¬ 
posed X-ray mission ATHENA [113] and GeV gamma- 
ray mission HERD [114] have made achieving similar 
energy resolution a priority, which will improve existing 
limits [115-126]. We encourage other missions to pursue 
this aggressively. 
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Supplemental Materials 


Outline: We first briefly discuss the mass models and 
dispersion profiles used to derive the results presented 
in the main text. We then provide an expanded discus¬ 
sion of two additional applications of DM velocity spec¬ 
troscopy, namely: probing the intrinsic DM dispersion 
profile using LOS observations and using velocity spec¬ 
troscopy of M31 to test detector causes. 

Radial Velocity Dispersion: To calculate the intrin¬ 
sic broadening of a DM line, a galactic mass model must 
be adopted to determine the velocity dispersion profile. 
Below, we describe the mass models used in our analysis 
of the MW and M31. 

Milky Way Mass Profile: We use model A1 of 
Ref. [47], which utilizes a DM halo determined by adia- 
batically contracting an initial NFW profile in the pres¬ 
ence of baryons. We summarize key aspects of the model. 

Before contraction, an NFW profile with scale radius 
Ts = 21.5 kpc is assumed to coexist with three axisym- 
metric baryonic profiles roughly associated with the nu¬ 
cleus, bulge/bar, and disc of the galaxy. The total bary¬ 
onic mass within a given radius is determined by the 
integration of the density profiles, with the addition of 
a central black hole of mass meH = 2.6 x 10 ®Mq. The 
enclosed baryonic mass is 

Mb(r) = mBH+/ [ dr'dnpi,{r')r'‘^. (8) 

Jo J 477 

The final DM profile is determined by contracting the 
initial NFW profile in the presence of this baryonic mass 
distribution. The baryonic profiles are adiabatically con¬ 
tracted under the assumption that spherical shells of 
matter do not cross and that the DM particles follow cir¬ 
cular orbits. This deepens the potential well and causes 
the DM to contract. Angular momentum conservation 
then dictates the following equations: 

G [Mb(r/) -I- Mdm(?■/)] rf = G Mhaio{?'i) n (9) 

' ^dm 

where Mhaio(?'i) is the halo mass before contraction and 
ridm and rib are the dark and baryonic matter densities, 
taken to be in the ratio ndm/(nb -I- Ddm) = 0.9; more 
recent observations give ndm/(nb -I- fldm) = 0.84 [86], 
which gives identical results. 

These equations are solved numerically to give a final 
radius, r/, corresponding to a given initial radius, r^. The 
contracted profile has a normalization p^ir = 8kpc) ~ 
0.4 GeV cm-3. 

The combined baryonic and contracted DM profiles are 
integrated to give the total mass enclosed within a given 
radius, Mtot(r): 

Mtot(r) = Mb(r) + Mdn>(r), (10) 


where M^^{r) is the dark matter mass within a radius r. 

The velocity dispersion is determined by the potential 
well of the galaxy, which is, in general, non-spherical. 
We approximate the true mass distribution by the spher¬ 
ically averaged mass profile given above. This approxi¬ 
mation has little impact outside of r ^ few kpc (where 
the DM becomes the dominant mass component), but 
greatly simplifies the calculation of the dispersion profile. 
Spherical symmetry allows for a simpler treatment of the 
Jeans equations [45] and, together with equilibrium and 
an isotropic velocity distribution, yields the expression 
for the radial velocity dispersion given in the main text. 

This mass profile (Eq. 10) generates a rotation curve, 
VciTcir) = y^GMtoti'f)/r, which is in good agreement 
with observations and a dispersion profile which agrees 
with results of previous papers [46, 47]. 

M31 Mass Profile: We use the mass model of 
Ref. [84]. Generalized Einasto profiles (given below) are 
used to describe the baryonic components 


Pb{a) = Pc exp 




( 11 ) 


with Pc, djy, Uc, and N adjusted to match data. The 
baryonic mass model include five components (nucleus, 
bulge, disc, young disc, and stellar halo). Together with 
the adopted NFW DM profile, the measured M31 rota¬ 
tion curve is reproduced well [84]. 

We also include a black hole of mass ttibh = 3.5 x 
lO^Mo [47]; more recent observations suggest a slightly 
larger mass or 1.4 x 10 ®Mq [82]. The inclusion of a cen¬ 
tral black hole yields larger velocity dispersions at small 
radii (< 10 pc), which increases the intrinsic width of DM 
lines arising from small angle LOS directions. However, 
because we focus on large angles {I ~ 10° — 40° in the 
MW and ~ 0.5° — 1.5° in M31), we do not probe the 
region affected by the black hole, so its effect is negligi¬ 
ble; we verified that our results were unmodified by this 
addition. Dispersions in M31 are comparable to those 
in the MW, but are systematically higher because of its 
larger mass and concentration. 

Figure Al shows the radial DM velocity dispersion pro¬ 
files for the MW and M31. Vertical bands represent the 
range of radii that contribute 90% to the signal along 
ip{l = 20°, 15] = 5°) in the MW and = 1° in M31. 

LOS Velocity Dispersion: The velocity distribution 
of DM is of great interest both for the information it con¬ 
tains about the particle nature of DM and for its implica¬ 
tions for direct and indirect detection experiments [55]. 
For example, models of self-interacting DM (SIDM) pre¬ 
dict higher velocity dispersions near the centers of DM 
halos. By measuring the LOS velocity dispersion, it may 
be possible to constrain SIDM interaction cross-sections, 







FIG. Al. Radial velocity dispersion profiles for the MW and 
M31. Shaded vertical bands indicate the range of radii that 
contribute 90% of the signal along = 20°,|fe| = 5°) in 
the MW and tp = 1° in M31; the radius ranges for the other 
directions discussed in the text are similar. Note that the 
lower bounds of these ranges are the smallest r probed by 
these directions. 

particularly in clusters where deviations between SIDM 
and CDM dispersions are large [76, 85]. Additionally, 
because sub-halos generate smaller velocity dispersions, 
variations in line width along different LOS could help to 
constrain the size and distribution of DM substructure. 

Because an observed DM signal will contain contri¬ 
butions from the entire LOS, and therefore a range of 
galactic radii, the full radial velocity dispersion cannot 
be probed directly. However, the observed LOS disper¬ 
sion may still contain useful information. It is natural 
to ask how well Astro-H may be able to reconstruct the 
intrinsic DM LOS dispersion, given the observed signal. 

Figure A2 shows both the intrinsic and observed (as¬ 
suming (TAH — 1-7 eV) LOS velocity width for a DM line 
in the MW and M31. Because the detector resolution is 
comparable to the intrinsic width, the detector response 
broadens the signal by a factor of ~ 

In principle, if the energy resolution of Astro-H were 
known exactly, the intrinsic width of the DM line could be 
reconstructed precisely; assuming the signal and detector 
response are both Gaussian, the effective width is a'^g = 
+ '^DM> can be determined simply. 

Of course, in practice, the resolution can never be 
known exactly. Assuming the goal uncertainty of 1 eV 
(the uncertainty is expected to be < 2 eV [71]), we esti¬ 
mate that the intrinsic width of a 3.5-keV DM line can 
be reconstructed with an uncertainty of ~ 40kms“^; for 



V’/'0s 

FIG. A2. Intrinsic (ctdm) and effective (nefif) LOS velocity dis¬ 
persion profiles for the MW and M31 as a function of -ipftps, 
the scaled angle relative to the center of each system. For 
the MW, ips = 50°, while for M31, tps = 2.5°; these scal¬ 
ings were chosen for display purposes. Intrinsic widths are 
determined by integrating the spectrum along the LOS using 
the radial velocity dispersion profiles given in the previous 
section. Effective widths include detector energy resolution. 
The increase in the LOS dispersion at small angles in M31 
is due to the rising radial dispersions shown in Fig. Al; for 
equally small (scaled) angles in the MW, the LOS dispersion 
decreases because only radii < 1 kpc, where the radial disper¬ 
sion is decreasing, contribute. 

higher energies the uncertainty in the width is smaller 
and scales as E~^. See the Appendix of Ref. [71] for 
more details regarding uncertainty in detector energy res¬ 
olution and intrinsic line width reconstruction. 

More speculatively, using information about the 
strength of the signal along the LOS, it may be possible to 
construct a course-grained radial velocity dispersion pro¬ 
file from the LOS dispersion. For example, we see from 
the vertical bands in Fig. Al that the range of radii that 
contributes to the ip{l = 20°,j6j = 5°) signal is narrow 
and that the dispersion of these points is directly reflected 
in the intrinsic LOS dispersion shown in Fig. A2. With 
additional pointings that probe different radii, it may be 
possible to constrain the radial dispersion profile using 
the measured line widths. This method would be most 
effective for small angles where the range of contributing 
radii is narrowest, although increased backgrounds would 
have to be overcome. 

Velocity Spectroscopy of M31: DM velocity spec¬ 
troscopy can also be applied to a signal observed from 
M31. Relative motion between the Sun and M31 pro- 
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FIG. A3. LOS velocity profiles for DM and HI gas [79] in 
M31. DM error bars are calculated assuming 2-Ms exposures 
with Astro-H and only CXB and detector backgronnds. 

duces a DM LOS velocity shift of ~ —300kms“^ that 
is essentially independent of viewing angle. For astro- 
physical lines, one must also consider the rotation of the 
M31 disc. This produces an additional LOS velocity shift 
that varies strongly with viewing angle, separating the 
DM and astrophysical lines by ± ~ 200kms“^ around 
±1° [79]. Detector lines are unshifted. 

The large differences in LOS velocities between DM, 
astrophysical, and detector lines make M31 a potentially 
powerful tool to probe the origin of spectral lines. How¬ 
ever, large LOS velocities are not by themselves sufficient 


to distinguish between these three causes; it is also neces¬ 
sary that the uncertainty in the profile centroid be small 
in comparison to the expected centroid separations. 

As discussed in the main text, the uncertainty in the 
centroid is given by SE = C{R) aes/V^t where tTeff 
is the observed line width, Ng is the number of signal 
events, and C(R) is a correction factor that accounts for 
the presence of backgrounds. As can be seen in Fig. A2, 
the observed widths of DM signals arising from M31 and 
the MW are expected to be quite similar. However, 
the number of DM signal events in M31 is considerably 
smaller, increasing the centroid uncertainty substantially. 

Figure A3 shows the LOS velocities for DM, astrophys¬ 
ical and detector lines as a function of the angular offset 
^ from the center of M31. We show the error bars on 
a DM signal assuming 2-Ms observations and only CXB 
and detector backgrounds. Astrophysical X-ray emission 
in M3I is not well studied outside of ~ 0.5°, but is ex¬ 
pected to be small [80, 81, 83]. However, even without 
including this background, it is clear that the significance 
(oc y/i) with which DM and astrophysical signals can be 
differentiated is considerably smaller than for the MW. 
However, the large differences between DM, astrophysi¬ 
cal and detector lines shifts could allow for cleaner sep¬ 
aration of these causes, if uncertainties were reduced. If 
MW observations of a line suggest a DM origin, several 
Ms would be well spent on M31 observations. 

Perhaps the greatest utility of observing M31 is in its 
power to test detector causes of a signal. This can be 
done most easily by looking directly at the center of 
M31. If the line is DM or astrophysical in nature, the 
signal strength should be strong and the centroid uncer¬ 
tainty correspondingly small, so that detector causes can 
be easily tested (Fig. A3). Though we have shown er¬ 
ror bars assuming 2-Ms observations, for this purpose, 
shorter exposures will clearly suffice. 





















